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The dipole-dipole interaction among ultra-cold Rydberg atoms is simulated. We examine a general 
interaction scheme in which two atoms excited to the x and x' states are converted to y and y states 
via a Forster resonance. The atoms are arranged in two spatially separated groups, each consisting 
of only one species of atom. We record the fraction of atoms excited to the y' state as the distance 
between the two groups is varied. With zero detuning a many-body effect that relies on always 
resonant interactions causes the interaction to have a finite range. When the detuning is greater 
than zero, another many-body effect causes a peak in the interaction when the two groups of atoms 
are some distance away from each other. To obtain these results it is necessary to include multiple 
atoms and solve the full many-body wave function. These simulation results are supported by recent 
experimental evidence. These many-body effects, combined with appropriate spatial arrangement 
of the atoms, could be useful in controlling the energy exchange among the atoms. 

PACS numbers: 32.80.Ee,37.10.Gh,03.67.Lx,02.70.-c 
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I. INTRODUCTION 

Mesoscopic ensembles of cold Rydberg atoms provide 
an ideal laboratory for exploring quantum dynamics. 
In the presence of a Forster resonance the dipole- 
dipole interaction allows for resonant energy exchange 
among the Rydberg atoms. Due to the large dipole mo- 
ments of Rydberg atoms, the interactions are long-range 
and take place on experimentally reasonable time scales. 
The interactions can be controlled in a number of ways, 
including manipulating the spatial arrangement of the 
atoms 0, [H, H[ , tuning a static electric field to shift the 
states of the atoms into resonance, and tailoring the mix- 
ture of Rydberg states [j| @. The potential for preci- 
sion control has led to a great deal of interest in using 
these systems for digital and analog quantum comput- 

ing 0,1! EM, Eld. 

In order to realize the potential of these systems, it is 
necessary to understand their complex many-body inter- 
actions. The energy exchange does not occur through a 
pairwise sum of binary interactions, but rather through 
the simultaneous interactions among many atoms. These 
many-body effects were first revealed through a broaden- 
ing of the resonant energy exchange, which could not be 
accounted for by simply considering two-body interac- 
tions [3, EH, EH] • Due to the potentially large number of 
atoms involved, simulation has proven to be a fruitful av- 
enue for understanding these systems [1, ES EM El) H3] • 
Previous work has shown that it is necessary to in- 
clude the full many-body wave function for as many as 
9 atoms in the calculations to reproduce experimental 
features 0, El, El H3 . 

We simulate the dipole-dipole interaction among Ryd- 
berg atoms for the four state system shown in Fig. QJa). 
By tuning the electric field, the states can be shifted by 



the Stark effect such that the energy gap E x — E y and the 
energy gap E y > — E x > are made equal. This field tuned 



Jy 

resonant interaction is 



x + x — ► y + y 



(1) 



where x —> y with dipole moment \i and x' — > y' with 
dipole moment v. There are also the always resonant 
interactions 



- y 
y' 
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(2) 



It has been suggested that the always resonant interac- 
tions contribute to an enhancement of the resonant inter- 
action [3, EH- A system of this type has been extensively 
studied in Rubidium [3, EH , where 

y = 24p 1/2 

x = 25si/ 2 
x = 33si/ 2 



V 



34p 3/2 



(3) 



In a recent report 0, van Ditzhuijzen et al. observe the 
spatially resolved dipole-dipole interaction between two 
groups of Rydberg atoms with the energy levels 
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49p 3/2 
49s 1/2 
41d 3/2 
42p 1/2 . 
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To facilitate comparison to experiment, many of the sim- 
ulations in this report are performed with parameters 
similar to the experimental values in Ref. (2j. 
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FIG. 1: (a) Energy level diagram for the resonant energy ex- 
change of the xx'yy' system, (b) The geometry studied con- 
sists of two groups of x and x' atoms separated by a distance 
d. 



The simulation is performed by diagonalization of the 
full dipole-dipole Hamiltonian matrix H . The Hamilto- 
nian in atomic units is given by 



tl — > P^.P.,./,,/ — 1 PxyPyx' 



2_^, PxyPx'y' t->3 



' R 3 



Px'y'Py'x' t->3 



m— n 



where m and n refer to individual atoms within each 
state and the sum is over all atoms in each state. The 
operators p a b take an individual atom from state a to 
state b, where a and b are the states of Fig. QJa). We 
ignore any orientation or spin effects and approximate 
the dipole-dipole interaction coupling by p 1 v/R 3 nn where 
R m n is the distance between the two atoms (this is simi- 
lar to the Hamiltonian given in Ref. [Hj]). The first term 
in Eq. ([5]) is the field tuned interaction and the next two 
terms are always resonant interactions. The final term 
gives the diagonal elements with a detuning, or energy 
defect, A = (E x + E x >) - (E y + E y >). While it has been 
found that dipole-dipole interactions can lead to conse- 
quential atomic motion [13, [H[ , we assume that we are in 
the regime of a "frozen gas." We therefore assume that 
the atoms are stationary on the time scales and densities 
studied. We also simplify the calculation by treating the 
dipole-dipole interaction as a process that occurs after 
the excitation of the atoms to Rydberg states, with no 
overlap in time. 

The possible states \4>i) of the system are enumerated 
in the xx'yy' basis. The initial state is assumed to be en- 
tirely composed of x and x' atoms and all excited states 
accessible from this initial state are included. For montc 
carlo simulations, the atoms are randomly placed in two 
groups consisting of exclusively x or x' atoms and sepa- 
rated by a distance d as shown in Fig.[ljb). The results 
are typically averaged over hundreds of runs. The fea- 
sibility of the monte carlo simulations is limited by how 
quickly the number of states scales with the total number 
of atoms. For this report, we simulated up to 16 total 
atoms; for the case of 8 x and 8 x' atoms this yields 12,870 



basis states. However, we found that the results compare 
well to experiment when including 12 total atoms. 



II. SIMULATION RESULTS 

We examine the general behavior of the energy ex- 
change between groups of Rydberg atoms by randomly 
placing the x and x' atoms in two spherical regions (see 
Fig. Gib)). The two spherical regions have the same ra- 
dius and are separated by a distance d. We simulate the 
interaction for some time (typically 10s of /is in 0.1 /xs 
steps) and the simulation results are averaged over 250 
runs at each value of d from /im (total overlap of the 
two regions) to 90 /jm in steps of 1 fim. We calculate the 
fraction of initial x' atoms found in the state y' , since this 
should be proportional to the experimental signal when 
using state selective field ionization. 

Fig. [2] shows the results for a detuning of zero when 
including 12 total atoms in the simulation. Six differ- 
ent cases were simulated: n x' atoms and (12 — n) x 
atoms where n = 1 ... 6. The dipole moments for each 
transition are equal, with fi = v = 1000 au. Due to 
this symmetry, the fraction of atoms in the y' state for 
n = 7 ... 1 1 can be inferred from the same data. For 
the case where the number of x and x' atoms are equal, 
shown in Fig.[2{a), the energy exchange persists to large 
separation. For the remaining cases, two of which are 
shown in Fig. [2Kb) and Fig. 0(c), the y' fraction drops to 
zero at smaller separation even though the detuning is 
zero. The finite range of the interaction in the absence 
of any detuning is due to a many-body effect that will be 
discussed in Sec. MI Al 

The most interesting feature appears when the detun- 
ing is greater than zero. Fig. [3J shows simulation data 
generated in the same way as the data in Fig. [2J with the 
exception that the detuning A w 2 MHz. In Fig. [3] the 
peak in the "strength" of the interaction, as measured by 
the y' fraction, is significantly away from the overlap of 
the two regions. The location of the peak is persistent 
for all times. In all cases shown in Fig. [31 the y' fraction 
drops to zero with a finite range due to the combined 
effect of the detuning and many-body effects. 

Before exploring the detailed behavior and origin of 
the peak in Sec. IIII CI we note that this feature has been 
seen by van Ditzhuijzen et al. in their recent work demon- 
strating interactions between spatially resolved groups of 
Rydberg atoms Q. We have run simulations roughly 
mimicking their experimental parameters; the results are 
shown in Fig. [5] The two regions of atoms are mod- 
eled as two gaussian beams each with a beamwaist of 
14 fj,m and a length of 250 /im. To account for differ- 
ent numbers of x and x' atoms in each beam, different 
cases for 12 total atoms that are similar to the Rydberg 
populations cited in the experiment are averaged. The 
detuning is A « 2 MHz. With some adjustment of the 
detuning, similar results can be obtained for numbers of 
atoms larger than 12. 
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10 20 30 40 50 60 70 80 90 
separation (p.m) 

FIG. 2: (color online) Fraction of initial x' atoms in the y' 
state as a function of the separation between two spherical 
groups of randomly placed atoms for four different times. The 
detuning is zero and the total number of atoms is 12. The 
times are: 1 /is (solid blue), 4 /is (dot-dashed red), 7 fis (dot- 
ted yellow) , and 10 fis (dashed green) . (a) When the number 
of x and x' atoms are equal (6 each), the interaction persists 
to large distances. The inset shows data out to 250 /tm show- 
ing that the y fraction eventually reaches zero as the Rabi 
period increases. For the case of (b) 3 x' atoms and 9 x atoms 
and the case of (c) 1 x' atom and 11 x atoms, the fraction of 
atoms in the y' state drops to zero with a shorter separation 
(around 50 /mi). 




10 20 30 40 50 60 70 80 90 
separation (urn) 

FIG. 3: (color online) Fraction of atoms in the y' state as 
a function of the separation between two spherical groups of 
randomly placed atoms after an interaction time of 5 /is. The 
detuning A ~ 4 MHz and the total number of atoms is 12. 
The most prominent feature is the peak in the interaction 
that occurs away from overlap for configurations with 3 to 6 
x atoms. 



the experimental data, particularly on the two aformen- 
tioned prominent features. First, the fraction of atoms 
excited to the y' state drops rapidly to zero from about 
40 to 50 /im. Second, the fraction of atoms excited to 
the y' state peaks at a separation between the beams of 
about 20 /im. However, our predicted mixing fraction 
is less than the observed mixing fraction. It is possi- 
ble that this is due to ignoring the temporal overlap of 
the Rydberg excitation process and the dipole-dipole in- 
teractions, which has been found to increase the mixing 
fraction fl9j ]. 



III. THEORY 

Three Atom Model: One x and Two x Atoms 
with Zero Detuning 



The data in Fig. |4] is graphed for positive and nega- 
tive separations, corresponding to the beam of x atoms 
being displaced to either side of the the x' beam. While 
the simulation is manifestly symmetric about the over- 
lap of the two beams (0 /im), we graph Fig. 0] in this 
manner to facilitate comparison to the experimental re- 
sults in Ref. 0]. Our simulation data agrees well with 



When there are unequal numbers of x and x' atoms 
in the two groups, the y' fraction drops to nearly zero 
within a finite range (see Figs. HJa) and[SJb)). Inter- 
ference from always resonant interactions among atoms 
in each group suppresses the field tuned interaction be- 
tween the two groups and leads to less population trans- 
fer to the y' state. In the simulation it is possible to 
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separation (|_im) 

FIG. 4: Fraction of atoms in the y state after an interaction 
time of 5 /zs as a function of the separation between the two 
groups of Rydberg atoms with volumes defined by excitation 
lasers with gaussian beam profiles. The detuning is A ~ 
2 MHz and the total number of atoms is 12. The location 
of the interaction peaks away from overlap at about 20 /im 
and the sharp turn-off of the interaction from about 40-50 /im 
compare well to the data in Ref. @]. 




separation (p.m) 

FIG. 5: Fraction of atoms in the y state as a function of the 
separation between two spherical groups of randomly placed 
atoms after an interaction time of 5 /j,s. In this case, there 
are 3 x' and 9 x atoms. The solid blue line is for A=0 and is 
from the same data as Fig.[2jb). The same simulation was run 
with the always resonant interactions removed and the result 
is shown with the dashed red line. While the y fraction is 
smaller, it also persists to large separations, implicating the 
always resonant interactions as the cause for the finite range 
of the interaction at zero detuning. 



remove the always resonant terms from the Hamiltonian 
in Eq. [5J effectively turning off the always resonant in- 
teractions. The result, when applied to the case of 3 x' 
atoms and 9 x atoms (Fig. [^b)), is shown in Fig. [5] At 
small separations, with the always resonant interactions 
active, the energy exchange is enhanced and more atoms 
are found in the y' state. However, at larger separations, 
the energy exchange is strongly suppressed. With no al- 
ways resonant interactions the y' fraction, while smaller, 
persists to large separations. 
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FIG. 6: (a) One x atom at a distance d from two x' atoms, 
which are separated by a distance R. We hold the distance 
R fixed and consider what happens as d is varied, (b) The 
maximum probability of the three atoms being in the one of 
the states |0 e i) or [ e 2 ) vs. the separation d in units of R. 



The effect of the always resonant interactions among 
groups of x and x' atoms can be examined analytically 
with a simple three atom model with zero detuning. An 
x atom is placed a distance d from each of two x' atoms 
that are separated by a distance R. This geometry is 
shown in Fig.[5fa), where the x atom is labeled 1 and the 
x' atoms are labeled 2 and 3. The field tuned interaction 
between the x atom and an x' atom is given by u = 
/iu/d 3 . The always resonant interaction between the two 
x' atoms is given by v — v 2 jR? '. For simplicity, in the 
following analysis, we choose fi = v. If fx ^ v , this will 
only change the distance scale of the predicted behavior. 
Numerical calculations show that the following results 
are insensitive to the exact placement of the atoms. 

In the zero interaction basis, we write the states of the 
atoms as 

\4> g ) = \x) 1 \x') 2 \x') 3 

\<Pei) = \y)x\y%\x% 

\4> e2 ) = \y) 1 \x') 2 \y') 3 (6) 

where the subscripts on the right hand side are the atom 
labels. The initial state \4> g ) is connected to the "excited" 
states \4>ei) and \<fi e2 ) via field tuned interactions between 
atom 1 and either atom 2 or atom 3. The two excited 
states are connected to each other via the always resonant 
interaction between atoms 2 and 3. 

We can write the time-dependent Schrodinger equation 

as 

ic g = uc e \ + uc e 2 
ic e i = uc g + vc e2 

ic e2 = uc g + VC e i (7) 

where c g is the amplitude for \4> g ), c e \ is the amplitude 
for |0 e i) , and c e2 is the amplitude for |<^ e 2)- The solutions 
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when the atoms are initially in state \<p g ) are 

9 2V8u 2 + v 2yy ' 

+ (V8u 2 -v 2 + v )e^ uA+vA -^ t ) 



C e l 



C e 2 



2%/8w 2 + w 2 



( e -i(V8« 2 +f a +«)* _ e ^(V8u- ! +v , ' ! -v)t^ 

(8) 



We examine the behavior of the system in terms of P\2 = 
1 — |c 5 | 2 , the probability of finding the atoms in one of 
the excited states. From Eq. [5J we find 

i(VSu 2 + v 2 t). (9) 



Au 2 4 v 2 

P\2 — 1 — „ o , ~ + 



4u 2 



8u 2 4 v 2 8u 2 4 v 2 

Note that by setting the always resonant interaction 
v = 0, we recover standard Rabi oscillations whose only 
dependence on the separation d is in the frequency. 

The case where the x atom is placed equidistant be- 
tween the x' atoms (d = R/2 or u = 8v) yields a solution 
nearly identical to standard Rabi oscillations except that 
the maximum probability P\ 2 max is slightly less than one. 
As we increase d and move the x atom farther from the 
pair of x' atoms, the maximum probability of the system 
being found in one of the excited states \4>ei) or |0 e2 ) de- 
creases rapidly, as shown by plotting Pi 2m ax m Fig.[^b). 

As d/R increases, the always resonant interaction in- 
creasingly suppresses the field tuned interaction, even at 
zero detuning. This effect is similar to the dark states 
phenomenon [HI]. The dressed states of this three atom 
system, which are superpositions of the states of Eq. [51 
depend on d. As d increases, one of the dressed states be- 
comes nearly identical to \4> g ). The initial state \<p g ) thus 
becomes dark, totally decoupled from the other states. 




FIG. 7: Two x atoms and two x' atoms arranged so that the 
separation between any xx' pair is always d. The pair of x 
atoms and the pair of x' atoms each comprise a "group" of 
atoms. The separation within the pairs is R and is held fixed, 
(a) The two groups of atoms at their smallest separation; all 
of the atoms are in the same plane. This view is along the 
axis of separation, (b) A view of the atoms from the "side." 
From this view, atoms 3 and 4 are on top of each other and 
the distance R between them is foreshortened. 



x' atom is d. As before, the field tuned interaction is 
given by u = [iv/d 3 and the always resonant interaction 
by v — [i 2 /d 3 . 

In the zero-interaction basis, we write the states of the 
atoms as 



I0 9 > = 




|z'>4 


I0el> = 


\y)i\x) 2 \y') 3 




\<t>ei) = 


\y)i\ x )2\ x '}3 


v')a 


|0e3> = 


k)ily) 2 l2/>3 




0e4) = 
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\4>es) = 


|y)ily) 2 |y% 





(10) 



where the subscripts on the right hand side are the atom 
labels. We can write the time-dependent Schrodinger 
equation as 



The evolution of the initial state to a dark state plays an 
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important role in explaining the peak in the interaction 
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as well. 
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B. Four Atom Model: Two x and Two x' Atoms 


id e 4 


= UCg 4 


2(u + 


v)c e i 4 


- vc e2 4 


- VCe3 + 


UCe5 


with Zero Detuning 


ic e5 


= AUCg 


+ UC e i 


4" UC e2 


+ uc e3 


4 UCei 





(11) 



The three atom model does not address the results of 
Fig. [UJa), which show that the range of the energy ex- 
change at zero detuning for equal numbers of x and x' 
atoms is limited only by the interaction time (see partic- 
ularly the inset in Fig. Ufa)). We construct the four atom 
model shown in Fig. [7ja)-(b) to examine the asymptotic 
behavior of the energy exchange at large separations for 
equal numbers of x and x 1 atoms. Numerical calculations 
show that the following results are insensitive to the par- 
ticular choice of geometry so we choose an arrangement 
that allows us to obtain an analytical solution. Both the 
two x atoms and the two x' atoms are separated by a 
distance R. The distance between any x atom and any 



where c g is the amplitude for the initial state \<p g ) and 
c e i is the amplitude for the excited state \<p e i)- 

If we keep only terms that are first-order in u, the 
solution for c g when the atoms are initially in state \4> g ) 
is 
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(12) 



so that the probability of the system being found in the 
state \4> g ) is 



cos 2 (4wi) 



9u 2 sin 2 (4ut) 
1024w 2 ' 



(13) 



At large separations, the only dressed states coupled to 
\4> g ) are equal superpositions of \<f> g ) and \<f>e&). The al- 
ways resonant interactions have a neglible effect and the 
energy exchange resembles simple Rabi oscillations be- 
tween the two groups of atoms. 



C. Interaction Peak away from Overlap 

A three or four atom model is insufficient to model the 
enhancement in the interaction away from overlap. This 
is evident if we vary the number of atoms included in the 
calculation. Fig. [5] shows simulation data for different to- 
tal numbers of atoms, from 6 to 16, with equal numbers 
of x and x' atoms in each case. This data is generated in 
the same fashion as the data for Figs. [Hand [31 with the 
exception of the 14 and 16 atom data. Since the num- 
ber of basis states is large for these two cases (3,432 and 
12,870), the 14 atom data was generated with reduced 
averaging (100 runs per 1 /im) and the 16 atom data 
was generated with reduced averaging and reduced res- 
olution (50 runs per 3 /im). For different total numbers 
of atoms the radius of the spherical regions is adjusted 
to hold the density constant. The location of the peak 
depends strongly on the number of atoms included in the 
simulation. 

With an amporphous sample of atoms, no peak in the 
interaction is observed in the 6 atom data in Fig. [5) We 
will see that under certain conditions a weak effect may 
be observed in the 6 atom case. However, no effect is 
observed for 5 atoms or fewer. The presence of the peak 
in the interaction away from overlap is therefore an in- 
trinsically many-body effect, requiring at least 6 atoms 
to manifest. 

In order to explore the many-body interactions in more 
detail, we examine the results for a simpler geometry of 
fixed instead of random atom positions. The two groups 
of atoms are arranged in a regularly spaced linear array, 
as shown in Fig. [9fa), with the atoms spaced by 10 fim. 
Using this model, we can study the dependence of the 
peak position on the detuning. Fig. [TO] is an intensity 
plot of the time-averaged y' fraction as a function of the 
separation between two lines of 5 atoms and the detun- 
ing. The detuning is varied from MHz to 10 MHz, 
corresponding to a few orders of magnitude smaller to a 
few orders of magnitude larger than the dipole-dipole in- 
teraction matrix elements, depending on the separation 
between the lines. 

As the detuning approaches zero, the peak in the in- 
teraction moves to larger separations and the interaction 
eventually becomes constant for all separations at A = 0. 
In an experiment, sources of detuning might include not 
only the experimental choice of electric field, but also 
the effects of electric field inhomogeneity and potentially 
a trapping magnetic field. In our simulations, we choose 
A to achieve reasonable agreement with experiment, as 
in Fig. [H in this sense A should be considered a free 
parameter. However, the location of the peak changes 
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separation (urn) 

FIG. 8: (color online) Fraction of atoms in the y' state as a 
function of the separation between two amorphous spherical 
groups of atoms. The detuning A m 4 MHz. Shown here are 
the cases for which the number of x and x' atoms are equal for 
different total numbers of atoms. The peak in the interaction 
away from overlap moves to larger separations and becomes 
smaller as the total number of atoms is increased. Unless a 
sufficient number of atoms is included in the calculation, the 
peak is not observed. 

rapidly only at small A. Once A as 2 MHz the peak lo- 
cation changes slowly, so our results are robust for a wide 
range of dctunings. 

The presence of the peak in the interaction away from 
overlap can be understood by examining the dressed 
states of our solution. In the non-interacting basis, we 
can write the states of the system as 

\<t>i) = \a) 1 \a) 2 ---\a) M \[3) M+1 \P) M+2 ---\P) N (14) 

where a = x or y, (3 = x' or y', the number of x and y 
atoms is M and the number of x' and y' atoms is (N—M), 
and the subscripts on the right hand side refer to indi- 
vidual atoms. The initial state \4>o) is composed of only 
x and x' atoms. At a particular separation d between the 
two lines, we can write the dressed states as a superpo- 
sition of the \4>i) with 

i 

The coefficient co,j determines the coupling of each 
dressed state to the initial state \4>o)- In Fig- M,c) the 
coefficients |coj| are plotted as a function of the sepa- 
ration between the two lines of atoms for the case of 5 
atoms in each line with a detuning of about 6 MHz and 
a spacing of 10 /im (see Fig.OJa)). As in the three atom 
model, one of the dressed states (solid blue) becomes 
nearly identical to the initial state beyond about 30 /im. 
The initial state is then dark, entirely decoupled from the 
other states (24|. As shown in Fig. [3(b), this corresponds 
to a region where the y' fraction rapidly approaches zero. 
However, at small separations (less than 20 /im), the 
eventual dark state is weakly coupled to the initial state. 
Another dressed state (dashed red) is strongly coupled 
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FIG. 9: (color online) (a) Geometry for a linear array of atoms 
where d is the distance between the lines and s is the spacing 
between the atoms, (b) Fraction of atoms in the y' state 
averaged over 10 j^s as a function of the separation d between 
two lines of 5 atoms with s=10 /im. The separation d starts 
at 1 /im and is increased in steps of 0.01 jj,m. As before, there 
is a peak in the interaction away from the overlap at about 
25 /»m. (c) The coefficients \co,j\ linking each dressed state 
to the initial state as a function of position. Only two of the 
252 dressed states are significantly coupled to the initial state; 
they are shown in solid blue and dashed red (the remaining 
250 coefficients are shown in dotted lines) . One dressed state 
(solid blue) is weakly coupled to the initial state at small 
separations but evolves to be nearly identical to the initial 
state at large separations. Another dressed state (dashed red) 
is strongly coupled to the initial state at small separations and 
evolves into a state that is nearly totally decoupled from the 
initial state at large separations. Where these two dressed 
states cross, they are both equally coupled to the initial state 
and coupled to many other excited states. At this crossing, 
there is a strong peak in the interaction shown in (b). 



to the initial state at small separations and thus weakly 
coupled to the excited states, corresponding to a region 
in Fig. EJb) where the y' fraction is small but non-zero. 
At separations between 20 and 30 /ini, these two dressed 
states switch roles and the coefficients |co,j| cross. In this 
region, both dressed states are equally strongly coupled 
to the initial state and well-mixed with many excited 
states. This mixing yields more population transfer to 
the y' state, creating a peak in the y' fraction away from 
overlap in Fig.[9^b). 

Figs. fTlTa)-(d) show a series of plots similar to Fig. [9j 
including the results for 4, 5, 6, and 8 atoms respec- 




40 60 
separation (p.m) 



100 



FIG. 10: Fraction of atoms in the y' state averaged over 10 /xs 
as a function of the separation and the detuning. The atoms 
are arranged as in Fig.[9ja) . Lighter shades of gray correspond 
to larger fractions of atoms in the y' state. The location of 
the peak in the interaction away from overlap can be seen as a 
light band that curves from the upper-left to the lower-right. 
The resolution of this data is 1 /im and about 0.07 MHz. 



tively. In the 4 and 5 atom cases, shown in Figs. [TlTa) 
and (b), there is no crossover point. The dressed state 
which evolves to become identical to the initial state is 
the most strongly coupled to the initial state even at 
small separations. Thus, there is no peak in the inter- 
action in Figs. [TlT a) or (b). When 6 or more atoms are 
included in the simulation, as in Figs. ITTT c) and (d), the 
coefficients cross and a peak away from overlap begins to 
form. 



IV. CONCLUSION 

Our results show, in agreement with other recent work, 
that in order to correctly model the collective interac- 
tions among Rydberg atoms it is necessary to calculate 
the full many-body wave function. In Refs. and [ll|, 
the authors conclude that at least 4 or 5 atoms must 
be included to accurately model their experiment and 
that summing over binary interactions is not sufficient. 
Here, we find good agreement using 12 atoms and re- 
quire a minimum of 6 atoms to observe the interaction 
peak away from overlap. However, we note that while 
the full many-body wave function is necessary to accu- 
rately model experiment, the results in Figs. O and [TJJ 
suggest that only two of the dressed states play a domi- 
nant role. Given the number of atoms and possible states 
involved in the simulation, this is indicative of the col- 
lective nature of the interactions. It also suggests that 
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some simplification of the analysis may be possible by 
considering the atoms collectively. 



separation (urn) 



separation (p.m) 



While it is clear that precisely positioning the Rydberg 
atoms will yield some degree of control over their interac- 
tions, our results show that that significant control may 
be possible even with amorphous samples. The peak in 
the interaction away from overlap offers a one such av- 
enue for experimental control. In Fig. [TO] there is a large 
region of parameter space where the interaction is nearly 
zero. Separating two groups of atoms and adjusting the 
detuning so that there is no interaction, one could place a 
system in this region. By slightly adjusting the detuning, 
one could then switch from no interaction to a relatively 
strong level of interaction. The interaction could be sim- 
ilarly controlled by changing the position of one or both 
groups of atoms. 



FIG. 11: (color online) Dressed state coefficients |cq,j| and 
y' fractions plotted together for the cases: (a) two lines of 2 
atoms, (b) one line of 2 atoms and one line of 3 atoms, (c) two 
lines of 3 atoms, and (d) two lines of 4 atoms. With 5 or fewer 
total atoms, as in (a) and (b), the y fraction peaks at overlap. 
With 6 total atoms, as in (c), a peak in the interaction away 
from overlap just begins to form near a 6 /im separation. 
As more atoms are included, as in (d) and Fig. [5] the peak 
becomes more pronounced and moves to larger separations. 



This material is based upon work supported by the 
National Science Foundation under Grant No. 0653544. 
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